Framework of the Alu Subfamily Evolution in the Platyrrhine Three-Family Clade of Cebidae, Callithrichidae, and Aotidae.

The history of Alu retroposons has been choreographed by the systematic accumulation of inherited diagnostic nucleotide substitutions to form discrete subfamilies, each having a distinct nucleotide consensus sequence. The oldest subfamily, AluJ, gave rise to AluS after the split between Strepsirrhini and what would become Catarrhini and Platyrrhini. The AluS lineage gave rise to AluY in catarrhines and to AluTa in platyrrhines. Platyrrhine Alu subfamilies Ta7, Ta10, and Ta15 were assigned names based on a standardized nomenclature. However, with the subsequent intensification of whole genome sequencing (WGS), large scale analyses to characterize Alu subfamilies using the program COSEG identified entire lineages of subfamilies simultaneously. The first platyrrhine genome with WGS, the common marmoset (Callithrix jacchus; [caljac3]), resulted in Alu subfamily names sf0 to sf94 in an arbitrary order. Although easily resolved by alignment of the consensus sequences, this naming convention can become increasingly confusing as more genomes are independently analyzed. In this study, we reported Alu subfamily characterization for the platyrrhine three-family clade of Cebidae, Callithrichidae, and Aotidae. We investigated one species/genome from each recognized family of Callithrichidae and Aotidae and of both subfamilies (Cebinae and Saimiriinae) of the family Cebidae. Furthermore, we constructed a comprehensive network of Alu subfamily evolution within the three-family clade of platyrrhines to provide a working framework for future research. Alu expansion in the three-family clade has been dominated by AluTa15 and its derivatives.

Analysis of Simian Endogenous Retrovirus (SERV) Full-Length Proviruses in Old World Monkey Genomes.

Simian endogenous retrovirus, SERV, is a successful germ line invader restricted to Old World monkey (OWM) species. (1) Background: The availability of high-quality primate genomes warrants a study of the characteristics, evolution, and distribution of SERV proviruses. (2) Methods: Cercopithecinae OWM genomes from public databases were queried for the presence of full-length SERV proviruses. A dataset of 81 Cer-SERV genomes was generated and analyzed. (3) Results: Full-length Cer-SERV proviruses were mainly found in terrestrial OWM, and less so in arboreal, forest- dwelling monkeys. Phylogenetic analysis confirmed the existence of two genotypes, Cer-SERV-1 and Cer-SERV-2, with Cer-SERV-1 showing evidence of recent germ-line expansions. Long Terminal Repeat (LTR) variation indicated that most proviruses were of a similar age and were estimated to be between <0.3 and 10 million years old. Integrations shared between species were relatively rare. Sequence analysis further showed extensive CpG methylation-associated mutations, variable Primer Binding Site (PBS) use with Cer-SERV-1 using PBSlys3 and Cer-SERV-2 using PBSlys1,2, and the recent gain of LTR motifs for transcription factors active during embryogenesis in Cer-SERV-1. (4) Conclusions: sequence analysis of 81 SERV proviruses from Cercopithecinae OWM genomes provides evidence for the adaptation of this retrovirus to germ line reproduction.

Evolution of the guanylate binding protein (GBP) genes: Emergence of GBP7 genes in primates and further acquisition of a unique GBP3 gene in simians.

Guanylate binding proteins (GBPs) are major players in the host immunity, providing defense against bacterial and viral invaders. Multigene families may suffer different processes of evolution. Gene families related to the immune system usually follow the birth-and-death evolution process, where duplicated genes can be deleted, gain new functions or become non-functional. We analyzed publicly available primate GBP sequences and their genomic organization and observed that GBP7 genes appear to have emerged from a duplication of GBP4 and seem to be only present in primates. Furthermore, GBP3 genes are only present in Simiiformes and probably originated from GBP1 genes. Finally, a duplication event occurred in the GBP6 in Tarsiiformes and became functional which might also explain the duplication of GBP6 in New World monkeys and Cercopithecidae. Taken together, this study provides new knowledge on the evolution of GBPs in primates and suggests that a revision of the GBPs nomenclature is necessary.

Primate phylogenomics uncovers multiple rapid radiations and ancient interspecific introgression.

Our understanding of the evolutionary history of primates is undergoing continual revision due to ongoing genome sequencing efforts. Bolstered by growing fossil evidence, these data have led to increased acceptance of once controversial hypotheses regarding phylogenetic relationships, hybridization and introgression, and the biogeographical history of primate groups. Among these findings is a pattern of recent introgression between species within all major primate groups examined to date, though little is known about introgression deeper in time. To address this and other phylogenetic questions, here, we present new reference genome assemblies for 3 Old World monkey (OWM) species: Colobus angolensis ssp. palliatus (the black and white colobus), Macaca nemestrina (southern pig-tailed macaque), and Mandrillus leucophaeus (the drill). We combine these data with 23 additional primate genomes to estimate both the species tree and individual gene trees using thousands of loci. While our species tree is largely consistent with previous phylogenetic hypotheses, the gene trees reveal high levels of genealogical discordance associated with multiple primate radiations. We use strongly asymmetric patterns of gene tree discordance around specific branches to identify multiple instances of introgression between ancestral primate lineages. In addition, we exploit recent fossil evidence to perform fossil-calibrated molecular dating analyses across the tree. Taken together, our genome-wide data help to resolve multiple contentious sets of relationships among primates, while also providing insight into the biological processes and technical artifacts that led to the disagreements in the first place.

Alu master copies serve as the drivers of differential SINE transposition in recent primate genomes.

Alu elements, averaging ~300bp in length, are a family of primate-specific short intersperse nuclear elements (SINEs) with more than one million copies and contributing to ~11% of primate genomes. Despite mostly being shared among primates, our recent study revealed highly differential recent Alu transposition among the genomes of primates from Hominidae and Cercopithecidae families. To understand the underlying mechanism, we analyzed six primate genomes and revealed species- and lineage-specific Alu profile exclusively defined by AluY composition. Among all Alus from the 6 genomes, we identified 5401 Alu master copies with 99% being from the AluY subfamily. The numbers of Alu master copies are positively correlated to the number of AluY elements in the genomes with the baboon genome having the largest number of most recent Alu master copies at high activities, while the crab-eating macaque genome having a low number of Alu master copies with low activity. Furthermore, the expression level of Alu master copies is positively correlated with their transposition activity. Our results support the concept that Alu transposition in primate genomes is driven by a small number of master copies, the number and relative activity of which contribute to the differential Alu transposition in recent primate genomes.

The draft genome of mandrill (Mandrillus sphinx): An Old World monkey.

Mandrill (Mandrillus sphinx) is a primate species, which belongs to the Old World monkey (Cercopithecidae) family. It is closely related to human, serving as a model for human health related research. However, the genetic studies on and genomic resources of mandrill are limited, especially in comparison to other primate species. Here we produced 284 Gb data, providing 96-fold coverage (considering the estimated genome size of 2.9 Gb), to construct a reference genome for the mandrill. The assembled draft genome was 2.79 Gb with contig N50 of 20.48 Kb and scaffold N50 of 3.56 Mb. We annotated the mandrill genome to find 43.83% repeat elements, as well as 21,906 protein-coding genes. The draft genome was of good quality with 98% gene annotation coverage by Benchmarking Universal Single-Copy Orthologs (BUSCO). Based on comparative genomic analyses of  the Major Histocompatibility Complex (MHC) of the immune system in mandrill and human, we found that 17 genes in the mandrill that have been associated with disease phenotypes in human such as Lung cancer, cranial volume and asthma, barbored amino acids changing mutations. Gene family analyses revealed expansion of several genes, and several genes associated with stress environmental adaptation and innate immunity responses exhibited signatures of positive selection. In summary, we established the first draft genome of  the mandrill of value for studies on evolution and human health.

Nomenclature report 2019: major histocompatibility complex genes and alleles of Great and Small Ape and Old and New World monkey species.

The major histocompatibility complex (MHC) is central to the innate and adaptive immune responses of jawed vertebrates. Characteristic of the MHC are high gene density, gene copy number variation, and allelic polymorphism. Because apes and monkeys are the closest living relatives of humans, the MHCs of these non-human primates (NHP) are studied in depth in the context of evolution, biomedicine, and conservation biology. The Immuno Polymorphism Database (IPD)-MHC NHP Database (IPD-MHC NHP), which curates MHC data of great and small apes, as well as Old and New World monkeys, has been upgraded. The curators of the database are responsible for providing official designations for newly discovered alleles. This nomenclature report updates the 2012 report, and summarizes important nomenclature issues and relevant novel features of the IPD-MHC NHP Database.

Cooperative evolution of two different TEs results in lineage-specific novel transcripts in the BLOC1S2 gene.

BACKGROUND: The BLOC1S2 gene encodes the multifunctional protein BLOS2, a shared subunit of two lysosomal trafficking complexes: i) biogenesis of lysosome-related organelles complex-1 and i) BLOC-1-related complex. In our previous study, we identified an intriguing unreported transcript of the BLOC1S2 gene that has a novel exon derived from two transposable elements (TEs), MIR and AluSp. To investigate the evolutionary footprint and molecular mechanism of action of this transcript, we performed PCR and RT-PCR experiments and sequencing analyses using genomic DNA and RNA samples from humans and various non-human primates. RESULTS: The results showed that the MIR element had integrated into the genome of our common ancestor, specifically in the BLOC1S2 gene region, before the radiation of all primate lineages and that the AluSp element had integrated into the genome of our common ancestor, fortunately in the middle of the MIR sequences, after the divergence of Old World monkeys and New World monkeys. The combined MIR and AluSp sequences provide a 3' splice site (AG) and 5' splice site (GT), respectively, and generate the Old World monkey-specific transcripts. Moreover, branch point sequences for the intron removal process are provided by the MIR and AluSp combination. CONCLUSIONS: We show for the first time that sequential integration into the same location and sequence divergence events of two different TEs generated lineage-specific transcripts through sequence collaboration during primate evolution.

Conservation of cell-intrinsic immune responses in diverse nonhuman primate species.

Differences in immune responses across species can contribute to the varying permissivity of species to the same viral pathogen. Understanding how our closest evolutionary relatives, nonhuman primates (NHPs), confront pathogens and how these responses have evolved over time could shed light on host range barriers, especially for zoonotic infections. Here, we analyzed cell-intrinsic immunity of primary cells from the broadest panel of NHP species interrogated to date, including humans, great apes, and Old and New World monkeys. Our analysis of their transcriptomes after poly(I:C) transfection revealed conservation in the functional consequences of their response. In mapping reads to either the human or the species-specific genomes, we observed that with the current state of NHP annotations, the percent of reads assigned to a genetic feature was largely similar regardless of the method. Together, these data provide a baseline for the cell-intrinsic responses elicited by a potent immune stimulus across multiple NHP donors, including endangered species, and serve as a resource for refining and furthering the existing annotations of NHP genomes.

Comparative Genomics Analysis Reveals High Levels of Differential Retrotransposition among Primates from the Hominidae and the Cercopithecidae Families.

Mobile elements (MEs), making ∼50% of primate genomes, are known to be responsible for generating inter- and intra-species genomic variations and play important roles in genome evolution and gene function. Using a bioinformatics comparative genomics approach, we performed analyses of species-specific MEs (SS-MEs) in eight primate genomes from the families of Hominidae and Cercopithecidae, focusing on retrotransposons. We identified a total of 230,855 SS-MEs, with which we performed normalization based on evolutionary distances, and we also analyzed the most recent SS-MEs in these genomes. Comparative analysis of SS-MEs reveals striking differences in ME transposition among these primate genomes. Interesting highlights of our results include: 1) the baboon genome has the highest number of SS-MEs with a strong bias for SINEs, while the crab-eating macaque genome has a sustained extremely low transposition for all ME classes, suggesting the existence of a genome-wide mechanism suppressing ME transposition; 2) while SS-SINEs represent the dominant class in general, the orangutan genome stands out by having SS-LINEs as the dominant class; 3) the human genome stands out among the eight genomes by having the largest number of recent highly active ME subfamilies, suggesting a greater impact of ME transposition on its recent evolution; and 4) at least 33% of the SS-MEs locate to genic regions, including protein coding regions, presenting significant potentials for impacting gene function. Our study, as the first of its kind, demonstrates that mobile elements evolve quite differently among these primates, suggesting differential ME transposition as an important mechanism in primate evolution.

Convergence of human and Old World monkey gut microbiomes demonstrates the importance of human ecology over phylogeny.

BACKGROUND: Comparative data from non-human primates provide insight into the processes that shaped the evolution of the human gut microbiome and highlight microbiome traits that differentiate humans from other primates. Here, in an effort to improve our understanding of the human microbiome, we compare gut microbiome composition and functional potential in 14 populations of humans from ten nations and 18 species of wild, non-human primates. RESULTS: Contrary to expectations from host phylogenetics, we find that human gut microbiome composition and functional potential are more similar to those of cercopithecines, a subfamily of Old World monkey, particularly baboons, than to those of African apes. Additionally, our data reveal more inter-individual variation in gut microbiome functional potential within the human species than across other primate species, suggesting that the human gut microbiome may exhibit more plasticity in response to environmental variation compared to that of other primates. CONCLUSIONS: Given similarities of ancestral human habitats and dietary strategies to those of baboons, these findings suggest that convergent ecologies shaped the gut microbiomes of both humans and cercopithecines, perhaps through environmental exposure to microbes, diet, and/or associated physiological adaptations. Increased inter-individual variation in the human microbiome may be associated with human dietary diversity or the ability of humans to inhabit novel environments. Overall, these findings show that diet, ecology, and physiological adaptations are more important than host-microbe co-diversification in shaping the human microbiome, providing a key foundation for comparative analyses of the role of the microbiome in human biology and health.

Complete nucleotide sequence characterization of DRB5 alleles reveals a homogeneous allele group that is distinct from other DRB genes.

Next Generation Sequencing allows for testing and typing of entire genes of the HLA region. A better and comprehensive sequence assessment can be achieved by the inclusion of full gene sequences of all the common alleles at a given locus. The common alleles of DRB5 are under-characterized with the full exon-intron sequence of two alleles available. In the present study the DRB5 genes from 18 subjects alleles were cloned and sequenced; haplotype analysis showed that 17 of them had a single copy of DRB5 and one consanguineous subject was homozygous at all HLA loci. Methodological approaches including robust and efficient long-range PCR amplification, molecular cloning, nucleotide sequencing and de novo sequence assembly were combined to characterize DRB5 alleles. DRB5 sequences covering from 5'UTR to the end of intron 5 were obtained for DRB5*01:01, 01:02 and 02:02; partial coverage including a segment spanning exon 2 to exon 6 was obtained for DRB5*01:03, 01:08N and 02:03. Phylogenetic analysis of the generated sequences showed that the DRB5 alleles group together and have distinctive differences with other DRB loci. Novel intron variants of DRB5*01:01:01, 01:02 and 02:02 were identified. The newly characterized DRB5 intron variants of each DRB5 allele were found in subjects harboring distinct associations with alleles of DRB1, B and/or ethnicity. The new information provided by this study provides reference sequences for HLA typing methodologies. Extending sequence coverage may lead to identify the disease susceptibility factors of DRB5 containing haplotypes while the unexpected intron variations may shed light on understanding of the evolution of the DRB region.

The Origin and Population History of the Endangered Golden Snub-Nosed Monkey (Rhinopithecus roxellana).

The origin and population history of the endangered golden snub-nosed monkey (Rhinopithecus roxellana) remain largely unavailable and/or controversial. We here integrate analyses of multiple genomic markers, including mitochondrial (mt) genomes, Y-chromosomes, and autosomes of 54 golden monkey individuals from all three geographic populations (SG, QL, and SNJ). Our results reveal contrasting population structures. Mt analyses suggest a division of golden monkeys into five lineages: one in SNJ, two in SG, and two in QL. One of the SG lineages (a mixed SG/QL lineage) is basal to all other lineages. In contrast, autosomal analyses place SNJ as the most basal lineage and identify one QL and three SG lineages. Notably, Y-chromosome analyses bear features similar to mt analyses in placing the SG/QL-mixed lineage as the first diverging lineage and dividing SG into two lineages, while resembling autosomal analyses in identifying one QL lineage. We further find bidirectional gene flow among all three populations at autosomal loci, while asymmetric gene flow is suggested at mt genomes and Y-chromosomes. We propose that different population structures and gene flow scenarios are the result of sex-linked differences in the dispersal pattern of R. roxellana. Moreover, our demographic simulation analyses support an origin hypothesis suggesting that the ancestral R. roxellana population was once widespread and then divided into SNJ and non-SNJ (SG and QL) populations. This differs from previous mt-based "mono-origin (SG is the source population)" and "multiorigin (SG is a fusion of QL and SNJ)" hypotheses. We provide a detailed and refined scenario for the origin and population history of this endangered primate species, which has a broader significance for Chinese biogeography. In addition, this study highlights the importance to investigate multiple genomic markers with different modes of inheritance to trace the complete evolutionary history of a species, especially for those exhibiting differential or mixed patterns of sex dispersal.

The evolution of S100A7 in primates: a model of concerted and birth-and-death evolution.

The human S100A7 resides in the epidermal differentiation complex (EDC) and has been described as a key effector of innate immunity. In humans, there are five S100A7 genes located in tandem-S100A7A, S100A7P1, S100AL2, S100A7, and S100AP2. The presence of several retroelements in the S100A7A/S100A7P1 and S100A7/S100A7P2 clusters suggests that these genes were originated from a duplication around ~ 35 million years ago, during or after the divergence of Platyrrhini and Catarrhini primates. To test this hypothesis, and taking advantage of the high number of genomic sequences available in the public databases, we retrieved S100A7 gene sequences of 12 primates belonging to the Cercopithecoidea and Hominoidea (Catarrhini species). Our results support the duplication theory, with at least one gene of each cluster being identified in both Cercopithecoidea and Hominoidea species. Moreover, given the presence of an ongoing gene conversion event between S100A7 and S100A7A, a high rate of mutation in S100A7L2 and the presence of pseudogenes, we proposed a model of concerted and birth-and-death evolution to explain the evolution of S100A7 gene family. Indeed, our results suggest that S100A7L2 most likely suffered a neofunctionalization in the Catarrhini group. Being S100A7 a major protein in innate defense, we believe that our findings could open new doors in the study of this gene family in immune system.

Functional characterization of enhancer evolution in the primate lineage.

BACKGROUND: Enhancers play an important role in morphological evolution and speciation by controlling the spatiotemporal expression of genes. Previous efforts to understand the evolution of enhancers in primates have typically studied many enhancers at low resolution, or single enhancers at high resolution. Although comparative genomic studies reveal large-scale turnover of enhancers, a specific understanding of the molecular steps by which mammalian or primate enhancers evolve remains elusive. RESULTS: We identified candidate hominoid-specific liver enhancers from H3K27ac ChIP-seq data. After locating orthologs in 11 primates spanning around 40 million years, we synthesized all orthologs as well as computational reconstructions of 9 ancestral sequences for 348 active tiles of 233 putative enhancers. We concurrently tested all sequences for regulatory activity with STARR-seq in HepG2 cells. We observe groups of enhancer tiles with coherent trajectories, most of which can be potentially explained by a single gain or loss-of-activity event per tile. We quantify the correlation between the number of mutations along a branch and the magnitude of change in functional activity. Finally, we identify 84 mutations that correlate with functional changes; these are enriched for cytosine deamination events within CpGs. CONCLUSIONS: We characterized the evolutionary-functional trajectories of hundreds of liver enhancers throughout the primate phylogeny. We observe subsets of regulatory sequences that appear to have gained or lost activity. We use these data to quantify the relationship between sequence and functional divergence, and to identify CpG deamination as a potentially important force in driving changes in enhancer activity during primate evolution.

The Targeted Sequencing of Alpha Satellite DNA in Cercopithecus pogonias Provides New Insight Into the Diversity and Dynamics of Centromeric Repeats in Old World Monkeys.

Alpha satellite is the major repeated DNA element of primate centromeres. Specific evolutionary mechanisms have led to a great diversity of sequence families with peculiar genomic organization and distribution, which have till now been studied mostly in great apes. Using high throughput sequencing of alpha satellite monomers obtained by enzymatic digestion followed by computational and cytogenetic analysis, we compare here the diversity and genomic distribution of alpha satellite DNA in two related Old World monkey species, Cercopithecus pogonias and Cercopithecus solatus, which are known to have diverged about 7 Ma. Two main families of monomers, called C1 and C2, are found in both species. A detailed analysis of our data sets revealed the existence of numerous subfamilies within the centromeric C1 family. Although the most abundant subfamily is conserved between both species, our fluorescence in situ hybridization (FISH) experiments clearly show that some subfamilies are specific for each species and that their distribution is restricted to a subset of chromosomes, thereby pointing to the existence of recurrent amplification/homogenization events. The pericentromeric C2 family is very abundant on the short arm of all acrocentric chromosomes in both species, pointing to specific mechanisms that lead to this distribution. Results obtained using two different restriction enzymes are fully consistent with a predominant monomeric organization of alpha satellite DNA that coexists with higher order organization patterns in the C. pogonias genome. Our study suggests a high dynamics of alpha satellite DNA in Cercopithecini, with recurrent apparition of new sequence variants and interchromosomal sequence transfer.

Nonhuman primate models of human viral infections.

Humans have a close phylogenetic relationship with nonhuman primates (NHPs) and share many physiological parallels, such as highly similar immune systems, with them. Importantly, NHPs can be infected with many human or related simian viruses. In many cases, viruses replicate in the same cell types as in humans, and infections are often associated with the same pathologies. In addition, many reagents that are used to study the human immune response cross-react with NHP molecules. As such, NHPs are often used as models to study viral vaccine efficacy and antiviral therapeutic safety and efficacy and to understand aspects of viral pathogenesis. With several emerging viral infections becoming epidemic, NHPs are proving to be a very beneficial benchmark for investigating human viral infections.

Superoxide dismutase 1 is positively selected to minimize protein aggregation in great apes.

Positive (adaptive) selection has recently been implied in human superoxide dismutase 1 (SOD1), a highly abundant antioxidant protein with energy signaling and antiaging functions, one of very few examples of direct selection on a human protein product (exon); the molecular drivers of this selection are unknown. We mapped 30 extant SOD1 sequences to the recently established mammalian species tree and inferred ancestors, key substitutions, and signatures of selection during the protein's evolution. We detected elevated substitution rates leading to great apes (Hominidae) at ~1 per 2 million years, significantly higher than in other primates and rodents, although these paradoxically generally evolve much faster. The high evolutionary rate was partly due to relaxation of some selection pressures and partly to distinct positive selection of SOD1 in great apes. We then show that higher stability and net charge and changes at the dimer interface were selectively introduced upon separation from old world monkeys and lesser apes (gibbons). Consequently, human, chimpanzee and gorilla SOD1s have a net charge of -6 at physiological pH, whereas the closely related gibbons and macaques have -3. These features consistently point towards selection against the malicious aggregation effects of elevated SOD1 levels in long-living great apes. The findings mirror the impact of human SOD1 mutations that reduce net charge and/or stability and cause ALS, a motor neuron disease characterized by oxidative stress and SOD1 aggregates and triggered by aging. Our study thus marks an example of direct selection for a particular chemical phenotype (high net charge and stability) in a single human protein with possible implications for the evolution of aging.

The Relationship between the (In-)Stability of NORs and Their Chromosomal Location: The Example of Cercopithecidae and a Short Review of Other Primates.

Amongst Cercopithecidae, the species of the Cercopithecini tribe underwent a very active chromosome evolution, principally by fissions, which increased their chromosome number up to 72. In contrast, all the species of Papionini have fairly similar karyotypes with 42 chromosomes. In animals, nucleolus organizer regions (NORs) are generally considered as instable structures, which frequently vary in size, number, and location at both infra- and interspecific levels. Although in Cercopithecinae the NORs, involved in breaks, exchanges, and translocations, behave like fragile sites in somatic cells, their number and location appear to be very stable between species. Fluorescence in situ hybridization of a 28S rDNA probe on metaphase chromosomes displayed a unique interstitial location in either an acrocentric pair (in 12 species of Cercopithecini) or a metacentric pair (in 6 species of Papionini). A non-exhaustive survey of literature data on NOR location in other primates shows that numerical variations of the NORs principally depend on their location: most multiple NORs are in terminal positions, while almost all unique NORs are in interstitial positions. We propose that this correlation is the consequence of the selection against gametic imbalances involving the chromosomal material distal to the NORs, which is effective when they are interstitially, but not terminally, located. Thus, the consequences of the interstitial NOR instability for reproduction are essentially limited to their size variations, as observed in Cercopithecidae.

Genetic Diversity and Structure among Isolated Populations of the Endangered Gees Golden Langur in Assam, India.

Gee's golden langur (Trachypithecus geei) is an endangered colobine primate, endemic to the semi-evergreen and mixed-deciduous forests of Indo-Bhutan border. During the last few decades, extensive fragmentation has caused severe population decline and local extinction of golden langur from several fragments. However, no studies are available on the impact of habitat fragmentation and the genetic diversity of golden langur in the fragmented habitats. The present study aimed to estimate the genetic diversity in the Indian population of golden langur. We sequenced and analyzed around 500 bases of the mitochondrial DNA (mtDNA) hypervariable region-I from 59 fecal samples of wild langur collected from nine forest fragments. Overall, genetic diversity was high (h = 0.934, π = 0.0244) and comparable with other colobines. Populations in smaller fragments showed lower nucleotide diversity compared to the larger forest fragments. The median-joining network of haplotypes revealed a genetic structure that corresponded with the geographical distribution. The Aie and Champabati Rivers were found to be a barrier to gene flow between golden langur populations. In addition, it also established that T. geei is monophyletic but revealed possible hybridization with capped langur, T. pileatus, in the wild. It is hoped that these findings would result in a more scientific approach towards managing the fragmented populations of this enigmatic species.